Thresholds for epidemic spreading in networks 
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We study the threshold of epidemic models in quenched networks with degree distribution given 
by a power-law. For the susceptible-infected-susceptible model the activity threshold A c vanishes in 
the large size limit on any network whose maximum degree k max diverges with the system size, at 
odds with heterogeneous mean-field (HMF) theory. The vanishing of the threshold has nothing to 
do with the scale-free nature of the network but stems instead from the largest hub in the system 
being active for any spreading rate A > l/\/k max and playing the role of a self-sustained source that 
spreads the infection to the rest of the system. The susceptible-infected-removed model displays 
instead agreement with HMF theory and a finite threshold for scale-rich networks. We conjecture 
that on quenched scale-rich networks the threshold of generic epidemic models is vanishing or finite 
depending on the presence or absence of a steady state. 
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The heterogeneous pattern of a network can have dra- 
matic effects on the behavior of dynamical processes run- 
ning on top of it [l[ , in particular when the distribution 
of the number k of contacts (the degree of an clement or 
vertex) exhibits long tails, as expressed by a power-law 
degree probability with the asymptotic form P(k) ~ fc~ 7 
Q. An example that has attracted a great interest due 
to its practical real-world implications is the modeling of 
epidemic spreading on contact networks Q . The simplest 
of these models is the SIS model Q , in which each vertex 
(individual) can be in one of two states, either suscepti- 
ble, or infected. Susceptibles become infected by contact 
with infected individuals, with a rate proportional to the 
number of infected contacts times a given spreading rate 
A. Infected individuals on the other hand become healthy 
again with a rate that can be set arbitrarily equal to 
unity. The model allows thus individuals to contract the 
infection time and again, leading, in the infinite network 
size limit, to a sustained infected steady state for val- 
ues of A larger than an epidemic threshold A c . In the 
SIR model [J], on the other hand, infected individuals 
recover (or die) and cannot change further their state. 
No steady state is now allowed, but a threshold still ex- 
ists above which the total number of infected individuals, 
starting from a very small infected seed, reaches a finite 
fraction of the network. The analysis of these and other 
models performed via a mean- field theory modified 
to take into account the heterogeneity of the network 
substrate 0, Q , led to the far-reaching conclusion that 
topological fluctuations, as measured by the second mo- 
ment of the degree distribution (fc 2 ), can have profound 
effects in many types of dynamics 0, Thus, for ex- 
ample, in the SIS model, the threshold takes the values, 
at the mean-field level A c = (k)/(k 2 ). For a long-tailed 
degree distribution with power law form, the second mo- 
ment diverges for 7 < 3, and one obtains the remarkable 
result of a vanishing epidemic threshold in the thermo- 



dynamic limit. These results have led to the widespread 
belief in the distinction between scale-free networks with 
7 < 3, where topology is highly relevant, and scale-rich 
networks with 7 > 3, where dynamical processes exhibit 
an essentially homogeneous mean-field behavior. 

In this Letter, building on some results previously re- 
ported, we present evidence that this belief is not correct 
for the SIS model on quenched networks (i.e. networks 
whose adjacency matrix is fixed in time) and that the 
scale-free nature of the contact pattern has no crucial 
effect on the value of the epidemic threshold. We inves- 
tigate the physical origin of this result, its validity for 
generic network structures and its consequences. On the 
other hand we show that for the SIR model the picture 
is different, a zero threshold occurring only in scale-free 
quenched networks. 

While heterogeneous mean-field (HMF) theory is ex- 
act on annealed networks (i.e. networks whose adjacency 
matrix in fixed only in average @), results beyond HMF 
theory for the SIS process on quenched networks (QN) 
have appeared in different contexts and with various lev- 
els of rigor. Already in 2003, Wang et al. Q argued that 
the epidemic threshold on an arbitrary undirected graph 
is set by the largest eigenvalue Ajv of the adjacency ma- 
trix, 

A c = A-\ (I) 

see also d, § . The relevance of Eq. (PJ becomes evident 
when it is complemented with the results of Chung et 
al. [HI, who calculated the largest eigenvalue of the ad- 
jacency matrix for a class of finite graphs with degrees 
distributed according to a power-law, obtaining 

" = U© %>VF cM N)' (2) 

where N is the network size, k c is network cut-off or 
degree of the most connected node (averaged over many 



network realizations [11|), and Cj are constants of order 
1. The cut-off k c is a growing function of the network 
size for uncorrelated scale-free networks, taking the value 
fc, ~ N 1 ' 2 for 7 < 3 and k c ~ AT 1 /(7-i) for 7 > 3 
12| . For 7 > 3 the ratio of the moments is finite and 
it is clear that the largest eigenvalue is governed by k c . 
Noticeably this remains true also for 5/2 < 7 < 3, since 
in that range (k 2 )/(k) ~ fe^~ 7 < v 7 ^- Only for 2 < 
7 < 5/2 the largest eigenvalue is set by the moments of 
the degree distribution. Combining Eqs. (TTJ) and ([2]), the 
behavior of the threshold for the SIS model in a power- 
law distributed network is, for sufficiently large size, 




7 > 5/2 
2 < 7 < 5/2 



(3) 



see also [13[. Since k c grows as a function of N for any 
7, the consequence of Eq. ([3]) is remarkable: In any 
uncorrelated quenched random network with power-law 
distributed connectivities, the epidemic threshold for SIS 
goes to zero as the system size goes to infinity. This has 
nothing to do with the scale-free nature of the degree 
distribution: It is always true as long as the cut-off k c 
diverges. Remarkably the threshold goes to zero also for 
Erdos-Renyi graphs (although logarithmically slow) , for 
which a formula similar to Eq. ^ exists |l4|. Differ- 
ent approaches [l3|, EH EH , have also pointed out that 
in the thermodynamic limit, the system is active for any 
A > 0. These results, however, have gone largely unno- 
ticed within the statistical physics community. 

A first issue raised by Eq. J3|) concerns the fact that, 
as any critical point, the epidemic threshold is well de- 
fined only in the thermodynamic limit. In a finite system, 
the dynamics is always doomed to fall into the healthy, 
absorbing state, even far above the threshold, due to 
stochastic fluctuations. The threshold for a finite net- 
work of size N must therefore be intended as the value 
separating the regime A < A c for which the epidemics 
decays exponentially fast (so that the expected survival 
time is of the order t ~ ln(iV)) from the regime A > A c 
where the survival time grows exponentially with TV to 
some power, r ~ e N , with a > 0. 

To investigate the validity of these results we have 
performed numerical simulations of the SIS model on 
quenched scale-rich networks with 7 = 4.5 and minimum 
de gr66 k-min — 3 , built using the uncorrelated configu- 
ration model [l7| . In order to compare results with the 
predictions in Eq. ^ one must take into account that the 
actual maximum degree k max in each network realization 
is a random variable, with average value (k max ) = k c . In 
particular, in the case 7 > 3, one can see [l8| that both 
the mean and the standard deviation of k max scale as 
k c ~ /V 1 /(7-i) j implying that k max always shows large 
fluctuations for different realizations of the degree se- 
quence. Therefore, we first consider networks in which 
k max has a fixed value, equal to the mean k c numerically 
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FIG. 1: Density of active sites for long times (restricted to 
surviving runs) in the SIS model on QN as a function of sys- 
tem size N, for 7 = 4.5 and different values of the parameter 
A. Notice that the straight line for A = 0.03 is due to the fact 
that no density smaller than 1/N can occur. 



estimated for the chosen system size N. In Fig. [T] we 
plot the density p a , calculated only for surviving runs, 
as a function of N for different values of A [l9| . Should 
the transition occur at a fixed value of A, p s would go 
to a constant for A > A c , decay exponentially for A < A c 
and as a power-law exactly at the transition. A com- 
pletely different behavior is observed: for all values of A, 
the curves are bent upward, indicating that the system 
is active for any A. This excludes the presence of a finite 
threshold for diverging N. While Eq. (fTJ) holds for SIS 
on any graph, Eq. (J2J) wa s instead obtained for a specific 
network model (intrinsically correlated for 7 < 3 and un- 
correlated for 7 > 3 [HI)- For generic topologies, it is 
simple to show [ll[ that \fk~ c is a lower bound for the 
largest eigenvalue of the adjacency matrix. This allows 
to conclude that, unless the degree distribution is strictly 
bounded from above, the threshold for SIS on any graph 
vanishes in the thermodynamic limit. 

How generic are these results? Prakash et al. Q 
have recently argued that Eq. (fTJ) is valid for all epi- 
demic processes, regardless of their particular micro- 
scopic details. To check this claim, we consider the SIR 
model. At the HMF level, the threshold takes the value 
Xf R = (k)/[(k 2 ) - (k)] 0, [H and is therefore finite 
for scale-rich networks with 7 > 3. From the analysis 
in Eq. on the other hand, it should be vanishingly 
small in the large network limit, according to Ref. [9(. 
We have checked this possibility by performing numerical 
simulations of the SIR model on networks with 7 = 4.5 
and different values of N, with fixed k mc 



— (kmax)- I n 

this case, the HMF estimated threshold takes the value 
\SiR ~ 0.31, independent of the network size, while the 
predictions from Eq. © are X^ IR ~ 0.0567, 0.0796, and 
0.1118 for the different network sizes considered. In Fig. [2] 
we report the final density of infected individuals R as a 
function of the spreading rate A, starting from a single 
randomly chosen infected node. The prediction of HMF 
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FIG. 2: Total number R of infected individuals in the SIR 
model on QN of different size N as a function of the spreading 
rate A. Networks have 7 = 4.5. 



FIG. 3: Number N s of active nodes in surviving runs for 
values of A smaller than the threshold predicted by HMF 
(X C (H M F) = 0.138 . . .) as a function of kmax, compared with 
the same quantity for star graphs. Networks have 7 = 3.5. 



theory seems in this case to be much more accurate than 
Eq. ([3]), contrary to the generic claim made in Ref. 
The threshold remains finite in the large N limit. 

To understand the different behavior of the two mod- 
els we look at the origin of the incorrect HMF prediction 
for SIS in QN. From a mathematical point of view, the 
HMF approach is equivalent Q to replacing the QN with 
given adjacency matrix by an annealed network with 
an averaged adjacency matrix, a^- 0,(2^. In the uncor- 
rected case this matrix reduces to Ojj = kikj/[N{k)], 
which has a unique nonzero eigenvalue A^r = (k 2 }/(k). 
Hence the annealed network approximation destroys the 
detailed structure of the eigenvalue spectrum of QN and 
preserves the correct largest eigenvalue only for 7 < 5/2 . 
This basic feature, and not (as suggested in Ref. [16() 
the disregard of dynamical correlations, is at root of the 
inaccuracy of the HMF approach. A more physical in- 
sight comes from the analysis of a star graph with one 
center connected to k max leaves of degree 1. In this 
case the largest eigenvalue of the adjacency matrix is 
A w = \/k max which implies A c = l/>/k max . The same 
result can be easily recovered by writing the rate equa- 
tions for the probability p max (p\ ) for the center (leaves) 
to be active, namely p max = -p max + (1 - Pmax)piXk rnax 
and p\ = —p\ + (1 — pijpmax^- Imposing the steady state 
condition one finds 



\ 2 K 



1 



(1 + Xk max )X 



Pi 



(1 + X)k rnax X 



(4) 



and hence the threshold condition above. The mes- 
sage of Eq. Q for a generic quenched random graph is 
strong: Independently from all the rest of the system, 
for A > 1/ 1 \Jk max the subgraph composed by the node 
with degree k max and its neighbors is in the active state. 
This core of activity provides a self-sustained source of 
infection that, since in the full graph the neighbors of 
the hub are not leaves, can transfer the activity to their 
other neighbors and spread in this way the epidemics to 
a finite fraction of vertices. This is confirmed by Fig. [31 



showing the number of actives nodes in surviving runs, 
N s , on a full network with 7 = 3.5 and on a star graph 
with the same k max 25|- For k max < 1/A 2 the values of 



N s in the full network and in the star graph are com- 
parable: both systems are subcritical and the subgraph 
centered around the node with degree k max is where ac- 
tivity lingers before disappearing. For k max > 1/A 2 the 
star graph becomes active and N s becomes asymptoti- 
cally proportional to k max . In the full network instead, 
the asymptotic behavior is N s ~ k2^ x ~ N, indicating 
that the active state is endemic: the hub spreads the ac- 
tivity to a finite fraction of the whole system. Reaching 
the the fully endemic state requires larger systems for 
small A, but nothing changes qualitatively for any A > 0. 

Understanding the behavior of SIS allows also to un- 
ravel why things go differently for SIR. In the former case, 
the possibility for hubs to be reinfected multiple times, 
which allows the presence of a steady state, boosts their 
impact on the dynamics. In the case of SIR, on the other 
hand, high-degree vertices can only be infected once and 
this strongly limits their role in the dynamics. Based on 
this observation, it is natural to conjecture that epidemic 
models allowing a steady state, such as SIS, will lead to a 
null threshold in any infinite QN, while all models with- 
out a steady-state will conform with HMF theory, with 
a finite threshold on scale-rich topologies. 

The strong effect of the hub in the dynamics raises fur- 
ther issues on the SIS model. While fixing the value of 
kmax to its ensemble average leads to results consistent 
with the presence of a non-zero threshold in finite sys- 
tems, as implied by Eq. ([3]), if this constraint is relaxed, 
kmax has large sample to sample fluctuations leading to 
nontrivial consequences. In Fig. @] we explore the effect 
of this variability by comparing simulations performed 
at fixed A and N, and different values of k max . The 
growth of the activity density for increasing k max in- 
dicates that the relation between the threshold (or the 
largest eigenvalue) and cut-off k c , Eq. can in fact 
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FIG. 4: Decay of the activity density for the SIS model in 
networks with 7 = 4.5, N = 10 6 , A = 0.1 and changing k mo .x- 



be refined, and be expressed in terms of the actual max- 
imum degree, A c = \j \J k max [18|. However, the large 
variations of k max among different realizations of the 



network with the same 7 and TV do not wash away as 
iV diverges and severely hinder the determination of the 
threshold in simulations with unrestricted k max . As men- 
tioned before, for 7 > 3 the standard deviation of k max 
increases as the average value (k max ) ~ jyt/h- 1 ) [ig| . 
and there is always a large sample to sample variabil- 
ity. Hence an unrestricted sampling at fixed A unwar- 
rantcdly averages networks with different thresholds and 
effective time scales, some subcritical and some supercrit- 
ical, making impossible even to determine the presence 
of a well-defined steady state. This fact is exemplified in 
Fig. [3J where we plot for comparison the activity density 
obtained averaging over networks with a freely varying 
kmax- For 7 < 3 the situation depends on the way the 
network is generated, and in particular on the way the 
upper bound of the degree distribution M = TV 1 /" grows 
[22j. If uj = 2 (uncorrelated configuration model [171 ]) or 
larger, the quantity (k 2 ) / (k) becomes sharply peaked as 
iV grows [22j. If instead oj = 1, as in the normal con- 
figuration model [23l| . the ratio (k 2 )/(k) (and hence the 
threshold) wildly changes from realization to realization, 
with relative fluctuations diverging as jV 2 ' 3 ^ 7 ^ 7 ^ 2 )^ 7-1 ) 
[22j] . Notice that in the intermediate region 5/2 < 7 < 3, 
the average value of \/k max is larger than (k 2 )/(k) but, 
since fluctuations of the latter diverge, for some network 
realizations the actual threshold A c is much smaller than 
the value predicted by Eq. Wc conclude that, unless 
7 < 3 and lo > 2, no average epidemic threshold can 
be properly defined from a numerical point of view for 
networks with unrestricted k max . 

In summary, we have studied how the threshold for 
models of epidemic spreading on quenched scale-rich net- 
works behaves as their size grows. The threshold for SIS 
model always vanishes in the thermodynamic limit, due 
to the role of hubs. This bears no relationship, at odds 
with the predictions of HMF theory, with the divergence 



of the second moment of the degree distribution, which 
is finite. For the SIR model instead the threshold van- 
ishes only for scale-free topologies (either quenched or 
annealed), in agreement with HMF theory. Wc conjec- 
ture that these different types of behavior are generic 
for systems possessing (or not) a steady state. While 
the result of a vanishing threshold for SIS is exact on 
quenched networks, it is however of limited interest from 
an epidemiological perspective. The interaction patterns 
over which real diseases spread generally vary over short 
time scales (24[, and are therefore better described by 
annealed topologies [(§, for which HMF theory works by 
definition, and the threshold is finite for 7 > 3. From a 
statistical physics point of view, instead, our results open 
a promising path towards a better understanding of the 
scope and limits of HMF theory as a theoretical tool to 
analyze dynamics on heterogeneous networks. 
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